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Abstract 

We study the effective potential of three-dimensional O(A^) models. In sta- 
tistical physics the effective potential represents the free-energy density as a 
function of the order parameter (Helmholtz free energy), and, therefore, it 
is related to the equation of state. In particular, we consider its small-field 
expansion in the symmetric (high-temperature) phase, whose coefficients are 
related to the zero-momentum 2j-point renormalized coupling constants g2j- 
For generic values of A^, we calculate g2j to three loops in the field-theoretic 
approach based on the e-expansion. The estimates of g2j, or equivalently of 
f2j = 92j/94~^, are obtained by a constrained analysis of the series that takes 
into account the exact results in one and zero dimensions. 
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expansion. 
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I. INTRODUCTION 



According to the universality hypothesis, most features of continuous phase transitions 
do not depend on the microscopic details of the systems, but only on their global properties 
such as the dimensionality and the symmetry of the order parameter (see e.g. Ref. 
The 0(iV)-symmetric universality classes describe many three-dimensional systems charac- 
terized by short-range interactions and an A^-component order parameter. We mention the 
liquid- vapour transition in classical fluids {N = 1), the A-transition in superfluid helium 
(A^ = 2), the critical properties of isotropic ferromagnetic materials (A^ = 3) and of long 
polymers (A^ — 0). The case A^ = 4 is relevant for high-energy physics: it should describe 
the critical behavior of finite-temperature QCD with two flavours at the chiral-symmetry 
restoring phase transition 0. Universality implies that critical exponents, as well as other 
universal quantities, are the same for all models belonging to the same 0(A^)-symmetric 
class. Thus, the universal results obtained for a representative of a given class, such as the 
0(A^)-symmetric (f)'^ Hamiltonian, can be used to predict the critical behavior of all systems 
in the same class. 

The effective potential (Helmholtz free energy) is related to the (Gibbs) free energy of 
the model. Indeed, if Ma = {(pa) is the magnetization and H the magnetic field, one defines 



where Z{H) is the partition function (the dependence on the temperature is understood). 
The global minimum of the effective potential determines the value of the order parameter 
which characterizes the phase of the model. In the high-temperature or symmetric phase 
the minimum is unique with M = 0, while in the low-temperature or broken phase J-'{M) 
presents a fiat region around the origin 0, i.e. J-'{M) is constant for \M\ < Mq where Mq 
is the magnetization at the coexistence curve. 

In the high-temperature phase the effective potential admits an expansion around M = 0: 



The coefficients can be expressed in terms of renormalization-group invariant quantities. 
We introduce a renormalized magnetization 

, i{t,H = (}fM{t,Hf 

^ = — ^Mnro) — • 

where t is the reduced temperature, x and ^ are respectively the magnetic susceptibility and 
the second-moment correlation length obtained from the two-point function of the order 
parameter 0, i.e. 



J^(M) = MH - ^\ogZ{H), 



(1) 



OO I 



(2) 



(4) 




1 Jdx x^G(x) 
2d Jdx G{x) 
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Then one may write 



= ^mV^ + j:rn''-^''-''^g2,v'^- (5) 



where 



and 



Here m = 1/^, g2j are functions of t only, and d is the space dimension. In field theory ip 
is the expectation value of the zero-momentum renormalized field. For t — the quantities 
g2j approach universal constants (which we indicate with the same symbol) that represent 
the zero-momentum 2j-point renormalized coupling constants. A simpler parametrization of 
the small-field expansion of the effective potential can be obtained by performing a further 
rescaling 

ip = — ^—2;, 6) 
which allows us to write the free energy as 

AT=—A{z), (7) 

^w4^^ + ^^^ + E^^.A (8) 

r2, = 4^ J > 3. (9) 

The function A{z) is related to the equation of state. Indeed one can show that z oc t~^M, 
and that the equation of state can be written in the form 

H oc t^^^. (10) 
oz 

The small-field expansion of the effective potential provides the starting point for the de- 
termination of approximate representations of the equation of state that are valid in the 
whole critical region. This requires an analytic continuation in the complex t-plane in order 
to reach the coexistence curve from the symmetric phase This can be achieved by 

using parametric representations 1^-0], which implement in a rather simple way the known 
analytic properties of the equation of state (Griffith's analyticity). This idea was success- 
fully applied to the Ising model, for which one can construct a systematic approximation 
scheme based on polynomial parametric representations and on a global stationarity 
condition 0, leading to an accurate determination of the critical equation of state and of 
the universal ratios of amplitudes that can be extracted from it In view of an ap- 

plication of such an approach to O(A^) modelsQ with > 1, we decided to improve the 



For A*" > 1, this approach is more difficult because of the presence of the Goldstone singularity at 
the coexistence curve, which must be somehow taken into account by the approximate parametric 
representations considered. 
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estimates of the coefficients r2j appearing in Eq. (P). It is worth mentioning that a better 
determination of the equation of state, and therefore of the universal ratios of amphtudes 
such as the ratio of the specific heat amphtudes A'^ /A~, is particularly important in the 
case N = 2, which describes the A-transition in ''He. A recent Space Shuttle esperiment |10 



made a very precise measurement of the heat capacity of liquid helium to within 2 nK from 
the A-transition, obtaining extremely accurate estimates of the exponent a and of the ratio 
A'^/A~. These results represent a challenge for theorists, who, until now, have not been 
able to compute universal quantities at the same level of accuracy (see e.g. Refs. [3,P^lHT3[ 
for recent theoretical estimates of a and A'^/A~). 

Several approaches can be employed to investigate the 0(A^)-vector models, such as 
field-theoretical methods starting from the 0^ formulation of the theory 



n 



(Fx 



(11) 



lattice techniques performing high- and low-temperature expansions, Monte Carlo simula- 
tions, etc... In order to study the small-field expansion of the effective potential, we consider 
the field-theoretic approach based on the e = A — d expansion W^. We extend the series of 



r2j to O(e^) (corresponding to a three- loop calculation) for generic values of N , and obtain 



new estimates from their analysis. Earlier estimates of r2j based on the e-expansion |T5 
were obtained using O(e^) series for generic values of and 0{e^) series for the Ising model, 
which were derived from the scaling equation of state known to O(e^) for generic values of 
[ig and to 0{e^) for the Ising model |3,|18[. 



Since the e-expansion is asymptotic, one needs to perform a resummation of the series 
in order to obtain reliable estimates. This can be efficiently done by exploiting its Borel 
summability and the knowledge of the large-order behavior [|l^]. Moreover, one may exploit 
exact results for low- dimensional models and perform constrained analyses of the e-series. 
The basic assumption is that the zero- momentum 2j-point renormalized couplings g2ji and 
therefore the ratios r2j, are analytic and quite smooth in the domain 4 > > (thus 
< e < 4). This can be verified in the large- A^ limit [^,0. One may then perform a 



polynomial interpolation among the values of d where the constants r2j are known, and then 
analyze the series of the difference. As we shall see, the analysis of the O(e^) series of r2j 
leads to a substantial improvement of the estimates of the first few r2j with respect to earlier 
results obtained from their O(e^) series |]T3. As a by-product of our analysis we also obtain 
new estimates for the ffist few r2j in the two-dimensional O(A^) models. 

The zero-momentum four-point coupling g = plays an important role in the field- 

m 



which provides an accurate de- 
In this approach, any universal 



theoretic perturbative expansion at fixed dimension 
scription of the critical region in the symmetric phase, 
quantity is obtained from a series in powers of g ( (^-expansion) , which is then resummed and 
evaluated at the fixed-point value of g, g* (see e.g. Refs. [|l|). Accurate estimates of g* have 
been obtained by calculating the zero of the Callan-Symanzik /3-function associated to g (see 
e.g. Refs. |l8|, |22|p3| , p!9|p^ ,p| ) . These results have been substantially con&med by computa- 
tions using different approaches, such as e-expansion [^, high-temperature expansion (see, 
e.g., Refs. @J2^ and references therein), Monte Carlo simulations p^^8|, etc... In this paper 
we reconsider the determination of g* from its e-expansion. In Ref. pO[ we calculated it to 
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O(e^), but unfortunately the series published there contains a numerical mistake^. For this 
reason, we report here the correct series and the results of the new analysis. We anticipate 
that the changes with respect to the estimates reported in Ref. I^O] are very small. 

The field-theoretic method based on the ^f-expansion at fixed dimension has been recently 
considered in the calculation of the zero- momentum couplings g2j with j > 2: for the Ising 
model five-loop series p|,§j3^ are available, while for generic values of A^, and gs have 
been determined to four and three loops respectively pT|-p^. The effective potential has also 
been studied by approximately solving the exact renormalization-group equations, providing 
some estimates of the coefficients g2j For the Ising model accurate estimates have been 
obtained from the analysis of lattice high-temperature expansions, see e.g. Refs. P, |35| , |36 . 
For > 1, which is the case considered in this paper, only the high-temperature of qq has 
been computed however, the results of its analysis are rather imprecise. 

The paper is organized as follows. In Sec. |I| we present the O(e^) series of that we 
have calculated. In Sec. |T| we give the results of the analyses of the series of r2j, which are 
then compared with the available estimates obtained in other approaches. The appendix is 
dedicated to the calculation of the three-loop integrals involved in the computation of the 
zero-momentum ra-point irreducible functions. 



II. EXPANSION OF r2j TO 0{e^). 

In the framework of the e-expansion we have calculated, to three loops, the one-particle 
irreducible correlation functions at zero momentum 

r2, = r(2/2....«,.,(o,...,o). (12) 

The number of diagrams one has to evaluate to compute r2j increases with increasing j. For 
example, the three-loop one-vertex irreducible diagrams necessary to compute Fg, Fg, Fio 
are 16, 36, 64 respectively. In our calculation we employed a symbolic manipulation package 
developed in Mathematica. It generates the diagrams using the algorithm described in 
Ref. , performs the necessary index contractions to determine the N dependence of each 



diagram, and compute the corresponding integral according to the procedure described in 
App. 0. 

The coefficients r2j of the expansion of A{z) in powers of z can be written in terms of 
F2j as 

, _ 92, _ (2j)! {N + 2y-' T2,n-' 



gj-i 2^3^-'j\Ul=2iN + 2i) F 
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From the three- loop expansion of F2j, one can derive the series of r2j to O(e^). In the 
following we report the series of the first few r2j, i.e. rg, rg and rio, that we will analyze in 
the following section. Writing 



The numerical expressions for the Feynman graphs appearing in Ref. [20| have been checked in 
Ref. |2£]. All numerical estimates are in agreement, except that of the constant H. The correct 
value is given in the Appendix. 
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r2j = Y.^2j,ie\ (14) 



i=l 

we find 



5(26 + N) 

"^''^ - 6(8 + AT) ^^^^ 
98 + 33iV + 4Ar2 (-8 + 7 N + N^) 

'^6,2 /o , aaS ^ o /o , An3 



+ A^)^ 3(8 + Ar)^ 

5 (17264 + 9968 N + 2574 A^^ + 319 + 7 N^) 



6(8 + AT)' 

5 A (-2176 - 172 AT + 152 AT^ + 9 N^) 20 (Qi + 7^5 A) (A^ - 1) 
^ 3 (8 + AT)^ ^ 3{S + Nf 

640 (A^ - 1) Q2 20 (682 + 49A^-2iV^) C(3) 
^ (sTiVf (8 + iV)^ 

35(80 + N) 

'''' ^ - 18(8 + AT) (1^) 

_ 35 (31904 + 7610 AT + 578 A^^ + 3 A^^) _ 5600 A (-8 + 7A^ + A^^) 
^^'^ ~ 54 (8 + A^)^ 27(8 + A^)^ 

35 (-259712 - 112232 A^ - 16204 A^^ _ 422 A^^ + 13 ^y^) 

18(8 + A^)^ 

35 A (105472 + 72528 A^ - 384 A^^ - 469 A^^) 2800(Qi + 7gA) (1 - A^) 
^ 81 (8 + AT)^ ^ 27(8 + A^)^ 

85120 (1 - N) Q2 70 (-29824 - 3010 A^ + 29 A^^) ((3) 
^ 9 (8 + Nf ^ 9 (8 + Nf 

35(242 + AT) 

"^^•^'^ - 3(8 + AT) ^^^^ 
35 (2083280 + 453428 AT + 28580 A^2 + 63 A^3) 162400 A (-8 + 7 A^ + iV^) 

TlO 2 Q 1 Q 

108(8 + Ar)^ 27(8 + Ar)3 

_ 35 (157284800 + 62464976 A^ + 8716080 A^^ + 388468 A^^ - 110 A^^ + 27 A^^) 



^8,3 



108 (8 + A^)^ 

140 A (-739808 - 822816 A^ - 35058 A^^ + 3359 A^=^) 
^ 81(8 + Ar)^ 

81200 (Qi + 7b A) {N - 1) 1433600 (A^ - 1) Q2 
27(8 + Ar)^ 9 (8 + A^)^ 

140 (463924 + 65932 A^ + 1585 A^^) ((3) 
^ 9 (8 + N)^ 



where 
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A = 1.171953619344729445... (18) 
Qi = -2.695258053506736953... 
Q2 = 0.400685634386531428... 

For = 1 the above expressions reproduce the O(e^) series that can be derived from the 
O(e^) equation of state calculated in Refs. [T^Jl^. We could also have computed the O(e^) 



series of r2j for some additional j > 5. Since, as we shall see, with increasing j, longer 
and longer series are necessary to obtain acceptable three-dimensional estimates from their 
analysis, we decided not to go beyond tiq. 

We now report the series of the fixed-point values of the zero-momentum four-point 



coupling g = (74, which corrects that given in Ref. which was plagued by a numerical 
mistake in one of the three-loop integrals. In the framework of the e-expansion, we found 
convenient to consider the rescaled coupling g, defined by 

1 (N + 8) f d\ 
g = ———S — 2--U, 19) 
^ 2(47r)^/2 3 V 2; ^ ' ^ ^ 

and expand it in powers of e. Due to the 0(e^^) factor multiplying g in the definition of g, 
the 0{e^) of g* corresponds to the O(e^) of g*. The e-expansion of g* to O(e^) is given by 

r{e) = T.9ke'', (20) 

k=0 



with 



go = 1 (21) 

3(14 + 3Ar) 



9i 

92 
93 



(8 + Ar)2 

1224 + 520A^ + 58A^2 _ 2Ar3 12(22 + 5A^)C(3) A(62 + 13iV) 

(8 + Ny (sTiVp 3(8 + iV)2 

341312 + 225312 + 57572 A^^ + 5494 iV^ - 99 ^ 4 A^^ (22 + 5 A^) vr^ 

8 (8 + A^)^ 15(8 + A^)^ 

2 (-3880 - 772 A^ + 431 A^^ + 90 A^3) ^(3) 43 (135 + 55 + 2 A^^^ ^(5) 

^ (sTiVf ^ (8 + N)* 
A(6500 + 2700Ar + 327N^ + 4N^) {Qi + 7^ A) (62 + 13 A^) 



2{8 + Ny Q{8 + Nf 

8(62 + 19A^)Q2 8H{22 + 5N) 



{8 + NY {8 + NY 



where 



H = -2.155952487340794361... (22) 
Finally we report the exact results in = 1 and d = that we will use in our constrained 



analyses of the e-series. In d = 1 we have [15 
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, 5NiN-in8N + 7) 

^ {N + 1){N + 4){4N - ly ^ ' 



_ 175 35A^(A^ - 1)2(256A^3 ^ 3037Ar2 ^ lygsAT - 588 
~ ^ 3(iV + l)(iV + 4)(A^ + 6)(4A^ - 1)3 ' 

rio = 1225 - 175iV(A^ - lf{N + 1)~-'^{N + 3)"^(iV + 4)-2(iV + 6)~^(A^ + 8)~^(4iV - 1)"^ x 



(149184 - 886968iV - 690826A^^ + 4219985A^^ + 6283975A^^ + 2913758A^' 
+552223A^^ + 44405A^^ + 1664A^''' 



for > 1, and 



re = 5, (24) 
175 



o 3 , 
no = 1225, 



for N < 1. In (i = we have 



_ 10(iV + 8) 

- 3(iV + 4)' 

_ 70(iV2 + 14A^ + 120) 



3(A^ + 4)(iV + 6) ' 
280(10752 + 3136iV + 256A^2 ^ 30iv3 + A^^^ 
(iV + 4)2(Ar + 6)(iV + 8) 



for N >1. It is not clear how to determine the value of r2j for = 0. Similarly to the d = 1 
case, one may conjecture that their values are independent of N for iV < 1, and therefore 
equal to those for = 1. 



III. ANALYSES OF THE SERIES. 

Since the e-expansion is asymptotic, it requires a resummation to get estimates for o? = 3, 
i.e. e = 1, which is usually performed assuming its Borel summability. The analysis of the 
series in powers of e can be performed by using the method proposed in Ref. |l^, which is 
based on the knowledge of the large-order behaviour of the series. It is indeed known that 
the n-th coefficient of the series behaves as ~ (— a)"r(n + 6o + 1) for large n. The constant a, 
which characterizes the singularity of the Borel transform, does not depend on the specific 
quantity; it is given by |]39| , ^ a = 3/(A^ + 8). The coefficient Bq depends instead on the 
quantity at hand. 

Consider a generic quantity R whose e-expansion is 

Rie) = Y.Rke'. (26) 

k=0 



According to Ref. [|1^], one generates new series Rp{a, b; e) according to 
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f poo 

R„{a, b]e) = y Bk{a, b) / dt e"* 



u{etf 
[1 - n{et)X 



(27) 



where 



u{x) 



a/1 + ax — 1 
\/l + ax + 1 



(28) 



The coefficients Bk{a,b) are determined by the requirement that the expansion in e of 
Rp{a, b; e) coincides with the original series. For each a, b and p an estimate of R is simply 
given by Rp{a, b;e = l). 

We follow Ref. |]20| to derive the estimates and their uncertainty. We determine an 
integer value of b, 6opt, such that 



-Rif (a, bopt] e = 1) ^ Rx^iia, &opt; e = 1), 



(29) 



for a < 1, where K is the highest order of the known terms in the series. In a somewhat 
arbitrary way we consider as our final estimate the average of RK{a,b;e = 1) with —1 < 
a < 1/2 and —2 + bopt ^ b < 2 + bopt- The error we report is the variance of the values of 
RK{a,b;e = 1) with — 1 < a < 1/2 and [&opt/3 — Ij <b< [46opt/3 + 1]. This procedure 
was already discussed and tested in Refs. |2^|T^. It seems to provide reasonable estimates 
and error bars, at least when a sufficiently large number of terms is known. Nonetheless, 
the method is ad hoc, and its reliability may depend on the quantity at hand. It is therefore 
important to check the accuracy of the error estimates for each quantity considered, for 
example by comparing, when it is possible, results obtained from different series representing 
the same quantity. 

It has been noted that, if one assumes a sufficiently large analytic domain in d for the 
quantity at hand, the estimates from its e-expansion can be improved using its known value 
in lower dimensions ETl-H3|,B0|,n!5|,f 



[]. One can check explicitly in the large- limit that the 

zero-momentum couplings g2j are analytic for < c? < 4. One may then conjecture that the 
analytic properties of g2j emerged in the large- iV limit are also valid for finite fixed values 
of A^. This point was already discussed in Refs. pil|JT3[ . 

Suppose that the values of the generic quantity R are known for a set of dimensions 
ei,...,efc. In this case one may use as zeroth order approximation the value for e = 1 of the 
polynomial interpolation through e = 0, ei,...,ek and then use the series in e to compute the 
deviations. More precisely, let us suppose that exact values i?ex(ei), • • •, -Rex(efc) are known 
for the set of dimensions ei, . . ., e^, > 2. Then define 



i=l 



e - e,; 



n 



(30) 



and 



(31) 



and finally 
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k 

i?imp(e) = [g(e) + 5(e)]n(e-Q). (32) 

i=l 

The resummation procedure is applied to S{e) and the final estimate is obtained by comput- 
ing -Rimp(e = !)• If the polynomial interpolation in d is a good approximation, one should 
find that the e-series which gives the deviations has smaller coefficients than the original 
one. Consequently one expects that also the errors in the resummation are reduced. In the 
cases considered, we find that, as expected, the coefficients of the series S{e) decrease in size 
with k, the number of exact values that are used to constrain the series. 

To begin with, we present the results of the analyses of the series for four-point coupling 
g* , cf. Eq. (pll). Table | shows the three-dimensional results for various constrained analyses, 
using the exact results m. d = 1 and c? = and the accurate estimates for two-dimensional 
models ||2^. Note that, for each the results of the various constrained analyses are 



consistent with each other. The two-dimensional estimates used in the constrained analyses 
are those already considered in Ref. [^], except for = 3. In this case we used the 
recent two-dimensional estimate g* = 1.7778(45) obtained from a form-factor bootstrap 
approach which should be more reliable than the high-temperature result g* = 1.724(9) 



considered in Ref. [^. For comparison, we note that for the two-dimensional 0(3) model 
our new estimate from the e-expansion is g* = 1.75(3). In Table |I| we present the two- 
dimensional results obtained by analyses constrained in d = 1 and d = 0. 

Tables H and |T| supersede the corresponding ones, i.e. Tables 1 and 2, of Ref. ||20|. The 



changes are small, so that the discussion presented there remains valid. In Table |Tl| we 
report the additional recent estimates which appeared after Ref. ||2^. These results have 
been obtained from: (i) a reanalysis of the six-loop /3-function in the framework of the 
fixed-dimension (yf-expansion (ii) the analyses of the high-temperature expansions (HT) 
of O(A^) cr-models on the cubic and bcc lattices |2^; (iii) high-temperature expansions of im- 
proved Hamiltonians (IHT) for which the leading scaling corrections are suppressed P, p^ , ^ |. 



The agreement among the various estimates is globally good. Let us only note that for = 3 
the error we obtained from our analysis of the e-series seems to be underestimated. A more 
complete list of references presenting estimates of g* can be found in Refs. p0|,|8|,|9|| . 

Since the series of r2j begins with a term of order e, we analyzed the O(e^) series of the 
quantity r2j/e. The constrained analyses were performed using the exact results of r2j for 
d = 1 and d = 0. We mention that in the case of the Ising model (A^ = 1) one may also 
use the precise two-dimensional estimates obtained from the analyses of the available high- 
temperature series for the Ising model on the square and triangular lattice These 
results were exploited in the analysis of the 0{e^) series of the Ising model presented in 



Ref. |T5|, allowing us to further improve the estimates of r2j. A discussion of the large- A^ 
limit of r2j can be found in Ref. |^ . 

Tables |IV| and |V| show respectively our three- and two-dimensional results for r^. There 
we report the estimates obtained from an unconstrained analysis and from analyses con- 
strained in various dimensions. Like the case of g*, the results of the various analyses are 
consistent with each other and the error decreases when additional lower dimensional values 
are used to constrain the analysis, supporting the estimate and the error obtained by the 
d = 0,1 constrained analyses. 

Tables |V| and |Vli] present respectively our three- and two-dimensional results for rg. In 
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this case, the results of the different analyses are not in complete agreement, indicating a 
possible underestimate of the uncertainty. Therefore, we believe that the final estimates of 
rg obtained from the d = 0,1 constrained analyses have an error larger than what is given by 
our algorithmic procedure. In Table |VIII| , where we report our final estimates, we multiply 
the errors by a factor of two: this quoted uncertainty should be more realistic. 

In the case of rio, the analyses, even those constrained, of its O(e^) series do not provide 
satisfactory estimates, but give only an order of magnitude. This may be explained by 
looking at the coefficients of the series of r2j: they increase rapidly with j, and, for rio, 
they are already much larger than its final three-dimensional estimate. We just mention 
that the d = 0,1 constrained analyses give rio = 29(34) for N = 2, rio = 16(24) for = 3, 
'"lo = 9(17) for = 4, where the reported errors are those obtained from our algorithmic 
procedure. 

Let us compare our results with the available estimates from other approaches. For 
7^ 1 there are not many published results: we are only aware of the estimates of Qq and 
gs presented in Refs. [^,^,0 (from which we can derive estimates of rg and rg using their 



results for g^). Table |VI1I| presents a summary of the available estimates of rg, rg and rio 
for several values of A^. For the sake of completeness, we report also the results for A^ = 1 
although the three-loop series were already known and no new results we have obtained in 
this work (a more complete list of results for the Ising model can be found in Ref. ||^). 
In Ref. pTI] g^ and gs were estimated from a Pade-Borel resummation of their d = 3 



(yf-expansion, calculated to four and three loops respectively. The authors of Ref. ||3T| argue 
that the uncertainty on their estimates of g^ is approximately 0.3%, while they consider their 
values for gs much less accurate. These results are in good agreement with ours, especially 
those for rg. We have also redone the analysis of the four-loop series of rg calculated in 
Ref. 1 511] using the same method employed here for the e-expansion (the results are reported 



in Table |VIII| ). In Table [VIII| we also show some results for the N = 2 model obtained from 



the analysis of the 20th-order high-temperature expansion of an improved lattice Hamilto- 
p6| : rg is in good agreement with the e- and (^-expansion, while for rg IHT seems to 



man 



give the most precise estimate. We also mention that Ref. |^| uses a renormalization-group 



approach (ERG) in which the exact RG equation is approximately solved (no estimates of 
the errors are presented there). Concerning rio, the result for N = 2 can be compared with 
the estimate coming from IHT: rio = ~13(7) P6[ . 

As already mentioned in the introduction, the estimates of the first few r2j that we have 
presented in this work may be very useful for the determination of the whole critical equation 



of state in 0{N) models with A^ > 1. Work is indeed in progress ^6|]. This approach was 
already successfully apphed to the Ising model. In the cases A^ > 1 the presence of the 
Goldstone singularity requires a more careful treatement. 
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APPENDIX A: THREE-LOOP DIAGRAMS. 



In the three-loop calculations of the zero-momentum irreducible 2j-point functions one 
deals with four families of integrals: 



B2{ni,n2,n3) = 

M(rai, 772,^3; 774, 77,5, Tie) = 



{2tiY (27r)^ 
d'^ki d<^k2 d'^ks 

d'^ki d% d% 



nz 



A{kir'A{k2rA{ki + k2) 



■A(fci)"^A(fc2)"^A(fc3)"^A(A;i + ^2 + ^3 



(Al) 



vn4 



{27tY {27tY {27lY 

X A(fci + A;2)"^A(A;i + k^ 
d'^ki d<^k2 d'^k^ 



A{k,rA{k2rA{k,) 



714 



(27r)^ (27r)^ (27r)'^ 
X A(A;2 - fc3)"^A(A;3)"^A(A;2)"« 

where A (A;) is the massive propagator, 

1 



A{k^)^'A{ki + k2r'A{ki + k3) 



A{k)^ 



P + 1' 



(A2) 



(A3) 



(A4) 



(A5) 



In order to compute these integrals, we have used an algorithm that expresses each quantity 
in terms of four basic integrals: -82(1, 1, 1), -83(1, 1, 1, 1), ^(1, 1, 1, 1, 1) and M(l, 1, 1, 1, 1, 1). 
The algorithm is based on the integration-by-parts technique [^9|-51|. The reduction of the 
integrals of type B3 is accomplished using the method presented in Ref. |5^, noting that 
these integrals are related to i?Ar(0, 0, 77i, 772, 71,3, 77,4), as defined in Ref. [|^. The integral 
53(1, 1, 1, 1) can be expressed in terms of the constant -B4 defined in since 



4 



B4 = — F 1 - 3e 

gSV J 3^4 



14 r(l - e/2)r(l + e)2r(l + 3e/2) 



r(H-e/2)2r(l + 2e) 



(1 e)(4 3e)(2 N.^B^ilAAAW - e/2)-'Til + e/2)-^ 



e2(2-e)2 



where d = 4 — e and 



(47r)^/2r(rf/2)' 

Using the expansion of B4 in powers of e reported in [^], we obtain 



(A6) 



(AT) 



53(1,1,1,1)= Nl 



5 1, 2^ 103 Stt^ ^ 



(A8) 



Let us now describe the algorithm for the other three cases. Let us begin with B2{ni, 772, 773). 
The basic relations are 
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{d - 2ni - ns)B2 = [-2ni 1+ - ^3 3+ + 3+(l" - 2'] 
ni-n2- n^)B2 = -{rii 1+ + ^2 2+ + 3+)52, 

I, n2, n^) = B2{ni ± 1, ?t,2, n^) and so on. Using (|A9 



{d — 111 — n2 — 



B2, 



(A9) 
(AlO) 



where l^i?2(?T-i, 'n.2, "^,3) = B2{ni ± 1, ?t,2, n3) and so on. Using (|A9|) and the relation obtained 
interchanging 2 and 3 we can reduce each B2{ni, n2, n^) to integrals of the form B2{mi, 1, 1) 
and to B2{mi,m2,0). The latter terms are the product of two one- loop integrals. To deal 



with the former ones, sum to (|A9D the relation obtained interchanging 2 and 3 in 
Then, use Eq. ( [A10| ) to eliminate the terms (^2 2+ + 3+)i?2 and {n2 2+ + 3"'")l~i?2- 
The new relation can be used to express each i?2(w^i, 1, 1) in terms of -82(1; 1? 1) and of the 
product of one-loop integrals. The integral B2{1,1,1) has been computed exactly in any 
dimension, obtaining p3| , ^ : 



52(1,1,1) 



2Nj 



e2(l-e)(2-e; 



r(l + e/2)2r(2-e/2)^ 



X k^'-'^'^^, + ^(1 - ^) 2Fr{l, e/2; 3/2; 1/4)" 



Expanding in powers of e we obtain 



^2(1, 1,1) 



..2(3 3 3 3A 
M 2e2 4e 4 4 8 



3 3,^ ^ ^ 



(All) 




where 



, 3-v/vr ^ v37r , , ^, 
A = -^Fi + ^(7i? + log3) 



6 
~[2 



(7i? + log3)^ 



r V3 



27r^ 



24 



(A13) 
(A14) 



il){x) is the logarithmic derivative of the F-function and Cl2(x) is Clausen's polylogarithm 
Here Fi and F2 are the following sums: 



F2 



1 ^ r(n + l)^(n+l) ^_„ 
2;to r(n + 3/2) 



4 V{n + 3/2) 



(A15) 
(A16) 



Let us now discuss the integrals (^(ni; ^2, n3; 714, ^5). First of all, note that if one of the 
indices is zero, G can be written as a B^ integral or as a product of a B2 and of a one- 
loop integral. We will now show that any G{ni\ n2, ris; n^, n^) can be expressed in terms of 
G'(l;l,l;l,l) and of G-integrals in which one of the indices is zero. Using the integration- 
by-parts technique, we obtain the following relations: 
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{d - 2ni - ng - 715)0 = \-{2ni 1+ + 3+ + 5^ 



+ n3 3+(r 



-(2^4 4^ 



{d — 2^4 — n5)G 

fSd 

n 

V2 

^5 — n4 + n4 4^ — n5 5^ 



2-) + 725 5+(l- -4- 
^5 5+) -ns 5+(l- - 



G, 



G, 



(A17) 
(A18) 



n2 - m - - n^j G = 
Hi 4+ (5" 



[ni 1+ + 722 2+ + ris 3+ + 724 4+ + 725 5+)G, (A19) 



5+(4" 



G= 0. 



(A20) 



First we consider Eq. (|A17|) and sum to it the three relations obtained interchanging (2^3), 
(23 ^ 45), and (23 ^ 54). Then, we use Eq. ( |A19|) ehminating (^2 2+ + ^3 3+ + 4+ + 
ns 5"'")G' and (n2 2"*" + ns 3"*" + 7x4 4+ + ng 5'^)1^G. Replacing rii by ni — 1, we obtain a 
relation which (for rii > 2) expresses G(ni; ^2, ns; n4, ns) in terms of G{mi; 1712, m^] 1714^, 1715) 
with nil < ^1- Therefore, in a finite number of steps we express the original integrals in 
terms of (7(1; n2, ^3; 214, ns) and of integrals in which one of the indices is zero. Then, we 
consider Eq. ( [A17| ) and use the relation obtained previously to eliminate l+G. The new 
relation can be used to reduce 77,3 to 1. Interchanging 2 and 3 we can similarly reduce n2 to 
1. We end up with G(l; 1, 1; 774,715) and with terms in which one index is zero. To further 
reduce the integrals we consider Eq. ( |A18| ). Replacing 714 by 714 — 1 we obtain a relation 
which reduces each integral to the form G{1; 1, 1; 1,715). Finally, we consider (A20) and use 
Eq. ( |A18|) to eliminate the terms 4'^G and 4~^5~G, obtaining a relation that reduces 715 to 
1. We should now compute G(l; 1, 1; 1, 1). Using the integration-by-parts technique we first 
obtain the exact relation 



^(1; 1, 1; 1, 2) = Go{l; 1, 1; 1, 2) - -Mo(l, 1, 1; 1, 1, 1). 



(A21) 



The integral Go(^i; ^^2, ""-s; ^^4, '^s) is given by Eq. ( |A^ ) replacing A(/c2)"^ with Ao(/c2)"^ 
and Mo(7ii, 7i2, 77,3; 774, 715, rie) by Eq. ( |A^ ) replacing A{k2 - ^3)"^ with Ao(fc2 - /^a)""*; Ao(A;) 
is the massless propagator l/k"^. For our purposes we only need the divergent part of 
Mo(l, 1, 1; 1, 1, 1) (the next contribution is reported in |]5B| ): 



Mod, 1,1; 1,1, 11 



2e 



[l + 0(e)]. 



(A22) 



The integral Go(l; 1, 1; 1, 2) can be computed exactly using the Mellin-Barnes technique 

(2-e)- 



57,58 . We obtain: 



Go(l;l,l;l,2) 



TC 



32sin=^(7re/2) (1 



+-(l-e)2Fi(l,6/2;3/2;l/4) + v^ ^ 



3^2(5/2, e,l;l-e/2, 1/2 + e/2; 1/4) 



2-l/2-e/2 2^-1 _ 



- 2' 



r (I) F (i±^) r (1 - 1) 



F(t)F(i + 
Using these expressions, we obtain finally for e 



Fi(e,3e/2;l/2 + e;l/4) 



(A23) 
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G(l;l,l;l,l) = iV3 




107]-^'? 7 1 

-y - y + ^(Qi + IeX + 2A) - SQ^ + ^C(3) + 0(e), (A24) 



where 



= ^ E I /O, — + 1/2) - ^(n + 1) + 2 log2 + - ) . (A25) 
4 r(?T, + 1/2) n! V ra/ 

Finally we consider M(ni, n2, na; n^, n^, ng). Again we will assume all indices to be positive: 
if one is zero, the integral can be rewritten as an integral of the G family. The basic relations 
we need are 

{d — 2nQ — Til — n^))M = 

i2ne 6+ + ni 1+ + ng 5+) + 5+(6" - 4") + rii l+(6" - 2")] M, (A26) 

+^5 5+(l" - 3") - ^3 3+(l" - 5") + ri4 4+(2" + 5" - 6" - 3")] M = . (A27) 



To compute this family of integrals we rewrite Eq. ( |A27 ) as 3+M = (^5 5^ + 



. . .)M. Then, replacing ^3 by ^3 — 1, we obtain a relation that reduces each integral 
to M(ni, n2, 1; ^4, ns, ng), or to M's in which one index is zero, and that, therefore, can 
be rewritten as G-integrals. Using the relations that are obtained replacing (123456) — *• 
(312645) and (123456 — > 321654) in Eq. ( [A27| ), we can similarly reduce n4 and 71,5 to 1. 
To further simplify the integrals, we consider Eq. ( |A27|) replacing (123456 231564) and 
rewrite it as M = l~{nQ 6"*" + . . .)M/ (ni — 1). If we apply this relation to M(ni, n2, 1; 1, 1, rig) 
we obtain three types of terms: (a) M(mi, 777-2, 1; 1, 2, 777,6) with 777i+m2+m6+l < 771+772+770; 

(b) M(777i, 7772, 1; 1, 1, tn%) with 7771 +"^2 + 7776 < 77i + 772 + (c) M (777i , 7772, 1; 1, 1, with 

777i + 7772 + = 77i + 772 + 776 and 777i < 77i. Tcrms of type (a) can be eliminated apply- 
ing repeatedly the relation we used to reduce 775, generating terms of type (b). Repeating 
the procedure, we end up with integrals in which either one index is zero or 771 is 1. An 
analogous procedure can be used to reduce 772 to 1. At the end of this reduction we should 
only consider M(l, 1, 1; 1, 1,776). We now consider Eq. ( |A26| ) and use the relations we con- 



sidered in the reduction of 775 and 771 to eliminate 5^M and 1^ M. Then, replacing ?7-g by 
77g — 1, we obtain M = 6^(. . .)M/{nQ — 1). This relation generates integrals of type G 
and l+6^6^Af, 5^6~6~M, and 6~M. The first two terms can be eliminated by applying 
repeatedly previous relations. At the end of the procedure any integral M is expressed in 
terms of M(l, 1, 1; 1, 1, 1) and of integrals of type G. The former integral has been computed 
in obtaining in d = 4 — e, 

M(l, 1, 1; 1, 1, 1) = Nl (^^C(3) + H + 0(e)) , (A28) 



where 



H . 4 Cli (0 . 2 L, (i) - ^ - g V 2 + i log^ 2, ( A29) 



15 



and Cl2(a;) and Li4(a:) are polylogarithms ||55|. A numerical determination of H was given 
in Ref. ||2^, rewriting 



H = ~ + jC(3)log5 + 3/?. 



(A30) 



20 



where H is the integral appearing in Eq. (B.19) of Ref, 
so that H —2.155953, in agreement with the numerical results of |]56[ 



Numerically H ^ -0.4346277, 
Unfortunately, in 



Ref. [^, we used H —0.04346277 — a factor of ten smaller — obtaining an incorrect 



estimate of H. 
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TABLES 



TABLE I. Three-dimensional estimates of g* from an unconstrained analysis, "unc" , and con- 
strained analyses in various dimensions. For the analyses which use the estimates in d = 2 we 
report two errors: the first one gives the uncertainty of the resummation of the series, the second 
one expresses the change in the estimate when the two-dimensional result varies within one error 



bar. 


N 


unc 


d=l 


d = 0,l 


d = 2 


d=l,2 


d = 0,l,2 





1.38(7) 


1.42(2) 




1.407(18+2) 


1.396(16+4) 




1 


1.40(8) 


1.44(3) 


1.41(3) 


1.424(22+0) 


1.408(19+1) 


1.408(12+1) 


2 


1.39(7) 


1.42(3) 


1.43(2) 


1.426(19+6) 


1.427(10+11) 


1.425(8+16) 


3 


1.39(7) 


1.40(3) 


1.41(2) 


1.410(19+1) 


1.420(7+2) 


1.426(7+2) 


4 


1.37(7) 


1.38(3) 


1.38(2) 


1.384(21+5) 


1.389(10+11) 


1.393(5+16) 


8 


1.31(5) 


1.30(3) 


1.30(2) 


1.301(19+1) 


1.304(9+1) 


1.307(4+2) 


16 


1.210(26) 


1.203(17) 


1.200(12) 


1.202(12+0) 


1.201(7+0) 


1.202(4+0) 


24 


1.160(17) 


1.155(12) 


1.151(9) 


1.152(9+0) 


1.150(5+0) 


1.150(4+0) 


32 


1.129(14) 


1.125(9) 


1.122(7) 


1.122(7+0) 


1.120(4+0) 


1.119(3+0) 


48 


1.091(10) 


1.089(6) 


1.087(4) 


1.087(4+0) 


1.085(3+0) 


1.085(2+0) 



TABLE II. Two-dimensional estimates of g* obtained from analyses constrained at d = 1 and 
at d = 0, 1. 



N 


d = l 


d = 0,l 





1.72(4) 




1 


1.84(6) 


1.76(5) 


2 


1.80(7) 


1.82(3) 


3 


1.73(8) 


1.75(3) 


4 


1.66(9) 


1.67(4) 


8 


1.46(6) 


1.46(3) 


16 


1.29(4) 


1.28(2) 


24 


1.21(3) 


1.20(2) 


32 


1.17(3) 


1.16(1) 


48 


1.12(2) 


1.11(1) 
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TABLE III. Three-dimensional estimates of g* = 5*(iV + 8)/(487r). The two results of Ref. fH] 
are relative to the cubic and bcc lattice respectively. 



N 


e-exp. 


d = 3 fif-exp. [§ 


HT H 


IHT 





1.396(20) 


1.413(6) 


1.388(5), 1.387(5) 




1 


1.408(13) 


1.411(4) 


1.408(7), 1.407(6) 


1.402(2) 


2 


1.425(24) 


1.403(3) 


1.411(8), 1.406(8) 


1.396(4) 


3 


1.426(9) 


1.391(4) 


1.409(10), 1.406(8) 




4 


1.393(21) 


1.377(5) 


1.392(10), 1.394(10) 





TABLE IV. Three-dimensional estimates of rg for various values of N from an unconstrained 
analysis of the e-expansion and constrained analyses m. d = 1 and d = 0, 1. The result in brackets 



for = 


has been obtained using 


our conjectured value for rg at d = 0, see Eq. ([25|). 




A^ 


unc 


d = 1 


d = 0,1 





2.180(80) 


2.148(22) 


[2.146(15)] 


1 


2.077(69) 


2.057(31) 


2.065(18) 


2 


1.980(65) 


1.955(28) 


1.969(12) 


3 


1.889(63) 


1.859(21) 


1.867(9) 


4 


1.812(66) 


1.778(23) 


1.780(8) 


8 


1.580(78) 


1.546(25) 


1.537(15) 


16 


1.333(38) 


1.310(17) 


1.300(18) 


32 


1.125(13) 


1.117(4) 


1.110(9) 


48 


1.036(10) 


1.033(2) 


1.029(4) 



TABLE V. 


Two-dimensional estimates of tq for 


various values of N from constrained analyses 


in d = 1 and d 


= 0, 1. The result in brackets for N 


= has been obtained using our conjectured 


value at d = 0. 






N 


d = 1 


d = 0,1 





3.745(47) 


[3.740(23)] 


1 


3.671(68) 


3.691(28) 


2 


3.494(58) 


3.530(18) 


3 


3.308(41) 


3.328(12) 


4 


3.155(44) 


3.159(12) 


8 


2.747(45) 


2.721(19) 


16 


2.368(34) 


2.335(24) 


32 


2.074(10) 


2.052(13) 


48 


1.950(4) 


1.937(6) 
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TABLE VI. Three-dimensional estimates of rg for various values of N from an unconstrained 
analysis of the e-cxpansion and constrained analyses in d = 1 and d = 0,1. The result in brackets 
for N = has been obtained using our conjectured value for rg at d = 0. 



N 


unc 


d= 1 


d = 0,1 





0.1(2.3) 


2.19(1.16) 


[3.13(53)] 


1 


-0.4(1.9) 


1.76(80) 


2.75(39) 


2 


-0.8(1.7) 


0.75(75) 


2.08(45) 


3 


-1.2(1.5) 


0.01(48) 


0.97(28) 


4 


-1.4(1.1) 


-0.50(31) 


0.19(23) 


8 


-1.8(4) 


-1.39(18) 


-1.18(9) 


16 


-1.7(4) 


-1.57(7) 


-1.54(4) 


32 


-1.3(1) 


-1.23(3) 


-1.24(2) 


48 


-0.97(5) 


-0.96(1) 


-0.969(4) 



TABLE VIL Two-dimensional estimates of rg for various values of N from constrained analyses 
in d = 1 and d = 0, 1. The result in brackets for A?^ = has been obtained using our conjectured 
value for at d = 0. 



A^ 


d=l 


d = 0,l 





24.6(2.3) 


[27.0(9)] 


1 


23.8(1.3) 


26.5(5) 


2 


19.7(1.4) 


23.2(6) 


3 


16.1(1.0) 


18.8(4) 


4 


13.5(8) 


15.4(3) 


8 


8.1(4) 


8.7(2) 


16 


5.0(2) 


5.1(1) 


32 


3.87(5) 


3.82(2) 


48 


3.71(2) 


3.64(1) 
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TABLE VIII. Three-dimensional estimates of rg and rg. When the original reference reports 
only estimates of g2j the errors we quote for r2j have been calculated by considering the estimates 
of g2j as uncorrelated. 



e-exp. d = 3 g-exp. HT ERG 



N 


= 


re 


2.148(22) [this work] 
2.1(3) H 


2.11(9) [this work] 












rs 


2(1) [this work] 

6(5) m 










N 


= 1 


re 


2.058(11) ll| 
2.12(12) li 


2.053(8) la 
2.060 iH 


2.048(5) i 
1.99(6) [H 
2.157(18) H 


1.92 1 


a 






rs 


2.48(28) H 
2.42(30) li 


2.47(25) la 
2.496 H 


2.28(8) i 
2.7(4) |35] 


2.18 


34] 






rio 


-20(15) PI 

-12.0(1.1) la 


-25(18) la 


-13(4) i 
-10(2) i 
-4(2) [H 






N 


= 2 


re 


1.969(12) [this work] 
1.94(11) II 


1.967 g 

1.970(40) [this work] 


1.951(14) H 
2.2(6) @ 


1.83 


34 1 






rs 


2.1(0.9) [this work] 

3.5(1.3) m 


1.641 H 


1.36(9) [H 


1.4 1 


1 


N 


= 3 


re 


1.867(9) [this work] 
1.84(9) (T| 


1.880 H 

1.884(32) [this work] 


2.1(6) 


1.74 1 


a 






rs 


1.0(0.6) [this work] 


0.975 




0.84 


34] 








2.1(1.0) m 










N 


= 4 


r& 


1.780(8) [this work] 

1.75(7) m 


1.803 [0 

1.809(27) [this work] 


1.9(6) 13 


1.65 


34] 






rs 


0.2(0.4) [this work] 

1.2(1.0) m 


0.456 




0.33 


g 
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